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1. Introduction 

Let C 5 C be a domain with boundary in an orientable smooth surface S 
in M^. Then the Gauss-Bonnet formula gives a method to compute the Euler-Poincare 
characteristic of D in terms of the Gauss curvature K of S, and the geodesic curvature 
kg of dD in S: 



2nx{D) = [ K + [ kg. 

J D JdD 



This can be found in any elementary book of Classical Differential Geometry. In this 
paper, we give a stereological version of this formula, in the following sense: let u E S'^ 
and denote by 7r„_A the plane orthogonal to u given by the equation {x, u) = A. When A 
varies on M, the different planes tTu^x can be considered as a "sweeping" plane in M.^ and 
its contact with D and with dD will give a method to compute x{D)'. 

Here, A is given by the points x E D such that D is tangent to Ttu,x at x for some A, B is 
equal to the subset of points y G dD where dD is tangent to tt^ ^ H S for some A, and k^ 
denotes the geodesic curvature of tt^ a H S. We will show that for a generic u E S^: 

1. the sets A, B are finite; 

2. if a; e A, then K{x) ^ 0; 

3. if 1/ e B, then 7r„ a fl S* is a regular curve in 5* and kgiy) ^ kg{y). 

In this way, the above formula has sense for almost any u E (that is, for any u E S'^\N, 
where is a null set in S"^). Moreover, this generalizes the result of 0, which computes 
the Euler-Poincare characteristic of a plane domain D C by looking at contact with 
lines. 



2. Proof of the formula 

The classical Poincare-Hopf Theorem stays that if S* is a closed orientable smooth 
surface and v is a smooth vector field on S with isolated zeros, then 

X{S) = J2 

v{x)=0 

where x{S) is the Euler-Poincare characteristic of S and the index, lndx{v), is just the 
local degree of v at x. As an immediate consequence, we get that if / : S* — > M is a Morse 
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function (that is, a function with non-degenerate critical points), then 

X{S) = ^ Ind.(/), 

xGS(/) 

where denotes the set of singular points of / and the index, lndx{f), is given by +1 
when X IS cL local extreme of /, or —1 when X IS cL saddle point. 

The Poincare-Hopf Theorem was generalized for surfaces with boundary by Morse in 
the following way: suppose that D is a compact orientable smooth surface with boundary 
(that we can assume embedded in an orientable smooth surface S without boundary) and 
let f be a smooth vector field on D such that 

1. V has isolated zeros on D, 

2. is not zero on dD, and 

3. V is tangent to dD only at a finite number of points. 

If 1/ G dD is a point where v is tangent to dD, generically we have that the integral curve 
of V at y is locally contained in S \ D or in D and we can assign an index +1 or — 1 
respectively. In the general case, we can extend this by using the Transversality Theorem 
and define an index, lndy{v) G {—1,0,+!}. Then, the Morse Theorem stays that 

xiD)= J2 Ind.(t;) + i 

v{x)=0 v{y)&TydD 

Now, we can interpret this in terms of critical points of functions. Suppose that / : 
D — i> M is a Morse function such that / has no critical points on dD and the restriction 
f\dD '■ dD ^ M is also a Morse function. For each point y G ^(/Ig/)), we have that the 
level set f^^{f{y)) is locally contained either in S* \ D or in D. In the first case, we 
say that y is a point of type "island" and we assign an index lndy{f) = +1 and in the 
second case, we say that it is point of type "bridge" and assign an index Indy(/) = — 1 
(see Figure 1). 




a) "Island" b) "Bridge" 

Figure 1 

Lemma 2.1. Let D be a compact orientable smooth surface with boundary and let f : 
D ^ be a Morse function such that f has no critical points on dD and the restriction 
f\dD '■ dD —>■ is also a Morse function. Then, 

X{D)= Ind,.(/) + i Yl 
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An interesting application of this formula is obtained when D G S CM.^ and we consider 
the family of height functions. Given u E S"^, we define the height function /z^ : — >■ M 
by hu{x) = {x, u). When we restrict this function to D (or S), the level sets are the plane 
curves 7r„ a H D (resp. 7r„ a fl S*), where 7r„^A is the plane h~^{\). In particular, if A is a 
regular value of hy\D (and hence of hy\s)i we have that vr^ ^ fl 5* is in fact a smooth curve 
in S. 

We will give a geometrical interpretation of the critical points oihu\D and hu\dD and the 
corresponding indices. Before that, we need to fix some notation. Given x G -D, we denote 
by N{x) and K{x) the normal vector and the Gauss curvature of -D at x respectively. 
Given y G dD, we denote by niy) and kgiij) the normal vector and the geodesic curvature 
of dD in D aX y respectively (we choose the orientation in dD so that n{y) points to the 
interior of D). Finally, we denote by k'^iji) the geodesic curvature of the curve vr^ ^ H S" at 
a regular point y of hu\D (in the case that vt^^a H S and dD are tangent at y, we choose 
in 7r„ A n S the same orientation). 

Theorem 2.2. Let D <Z S d he a domain with boundary in an orientahle smooth 
surface S in M'^ and let u E S^. 

1. Let X (z D. It is a critical point of huln if o-nd only if tIu,x is tangent to D at x 
(where hu{x) = \) if and only if u = ±N{x). 

2. Let X E D be a critical point of huln- It is non-degenerate if and only if K{x) ^ 0. 
Moreover, it is a local extreme of huln when K{x) > and a saddle point when 
K{x) > 0. 

3. Let y G dD be a regular point of hulc- It is a critical point of hu\dD if and only 
if T^u,\ r\ S is tangent to dD at y if and only if u = ±n{y), where it denotes the 
normalized orthogonal projection of u in TyD. 

4. Let y G dD be a regular point of huln which is also a critical point ofhy\dD- It is non- 
degenerate if and only if kg{y) ^ kg{y). Moreover, it is an island when kg{y) > k^{y) 
and a bridge when kg{y) < kg{y). 

Proof. Parts (1) and (2) are standard results on Generic Geometry and can be found for 
instance in [0]. We will prove (3) and (4). Suppose that a{s) is a parametrization of dD 
by arc length, so that y = a{so) and n{y) = N{y) A a'(so). Since {N{y),a'{sQ),n{y)} is 
a orthonormal basis of , we have that 

u = {u,N{y))N{y) + {u, a'{so))a'{so) + {u,n{y))n{y). 

In particular, provided that u ^ ±N{y), we have that u = ||, where 

V = {u,a'{so))a'{so) + {u,n{y))n{y). 

Thus, u = ±n{y) if and only if {u, a'i^So)) = if and only if y is a critical point of hu\dD- 
On the other hand, let (3{s) be a parametrization of 7^u,\ H 5 by arc length, with 
P{so) = y. In this case, {u,P{s)) = A for all s and thus {u,j3'{so)) = 0. This gives that u 
is also equal to the normal vector to P'{so) in TydD. In particular, y is a critical point of 
hu\dD if and only if the two curves are tangent. 

Suppose now that ?/ is a critical point of hu\dD- We have that 

a"{so) = kn{y)N{y) + kg{y)n{y), 

where kn{y) is the normal curvature of a. This gives that 

(M,a"(so)) = kniy){u,N{y)) + kg{y){u,n{y)). 
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But we can do the same with l3{s), getting 

= {u,P"{so)) = K{y){u,N{y)) + k^{y){u,n{y)) . 

This imphes that 

{u,a"{so)) = {kg{y) - kg{y)){u,n{y)) . 

Since {u,n{y)) ^ 0, we conclude that ?/ is a non-degenerate critical point of hu\dD if and 
only if kg{y) ^ kg{y). The fact that the cases kg{y) > kg{y) and kg{y) < kg{y) correspond 
to an island or a bridge respectively can be deduced from the way we have chosen the 
orientation on dD. □ 

The following consequence can be seen as a stereological version of the Gauss-Bonnet 
formula. 

Corollary 2.3. Let D G S G M.^ be a domain with boundary in an orientable smooth 
surface S in M^. For a generic u G 5^ the height function /i„|_D is a Morse function, has 
no critical points on dD and its restriction to dD is also a Morse function. In particular, 
we have that 

xiD)= Yl signK{x) + ^ ^ sign{k g{y) - k^{y)). 

xeD/N{x)=±u yedD/n{y)=±u 

Proof. The formula for x(-D) is a direct consequence of Lemma 3.1 and Theorem 3.2. We 
will see the first part of the statement. Given u G 5*^, h^lu is a Morse function if and only 
if u is a regular value of the Gauss map N : D S"^ (just note that the points x G D 
where K{x) = are the critical points of the Gauss map). Therefore, the Sard Theorem 
implies that for a generic u G S^, HuId is a Morse function. The other two conditions are 
also an easy consequence of this theorem. □ 

It is also interesting to look at some particular cases. Suppose for instance, that 5 is a 
plane, namely 5* = R^. Then we can use height functions : — > R for u G S^, instead 
of the height functions of M^. The level curves in this case are the orthogonal lines to u, 
which have zero curvature. Moreover, has no critical points in and the geodesic 
curvature is equal to the classical curvature of a plane curve. Thus, we get the following 
formula method to obtain the Euler number of a domain in based on 1-dimensional 
observations, modified from [^, and that has been used in Stereology, ||^. 

Corollary 2.4. Let D G M."^ be a domain with boundary. For a generic u G the 
restriction of the height function hu\dD is a Morse function and 

= ^ signA;(?/). 

y^dD /n{y)=±u 

In the case that 3 = 3"^, the level curves associated to u G S*^ are the corresponding 
parallels. Note that huls^ is a Morse function with two critical points, u and —u. Moreover, 
if we denote by 9u, Ju the spherical coordinates associated to u, we have that the geodesic 
curvature of the parallel 7„ = 70 is — tan 70. 

Corollary 2.5. Let D G 3"^ be a domain with boundary. For a generic u G 3'^ , 
^(^) = ^ H sign(A;3(|/)+tan7„(|/)) + #({M,-M}nZ)), 

y&dD /n{y)=±u 

where i^{{u, —u} fl D) means the number of times that u or —u belong to D. 



A STEREOLOGICAL VERSION OF THE GAUSS-BONNET FORMULA 



5 



It is also possible to compute the Euler number of D G from observations on a 
meridian that sweeps through D. This is not a consequence of the construction associated 
to the height functions, but a different family of functions. Let m G S*^ and let 5*^ be the 
circle in S'^ orthogonal to u. Then, we have the regular function 9^ '■ \ {u, —u} S^, 
whose level sets are the meridians 6^ = Oq (which have zero geodesic curvature). Moreover, 
for a generic u G 5*^, the restriction of 6u to dD is a Morse function, 0. 

Corollary 2.6. Let D G be a domain with boundary. For a generic u G S"^ , 

y£dD /n{y)A-U 

Proof. Let S = S"^ \ {u, —u} and let D be the domain with boundary in S given by 
D \ {Bi U B2), where Bi, B2 are small geodesic balls centered at u,—u which do not 
intersect dD. It follows from Lemma 3.1 that 

xiD) = ^ 5Z signkgiy). 
On the other hand, x{D) = x{D) + —u} Pi D), which concludes the proof. □ 

Remark 2.7. The preceding result may be extended to ovaloids in the following way: Let 
S be an ovaloid (compact and connected surface in for which the Gaussian curvature 
K > 0) and D a domain in S; then, S is diffeomorphic to a sphere S"^ through its Gauss 
map. Let D G denote the image of D under this diffeomorphism; then, x(-D) = x(-D) 
and it is possible to obtain x(-D) from observations in the curves which are image of the 
meridians in S"^. 

3. STEREOLOGICAL APPLICATIONS 

The Euler- Poincare characteristic for domains in and has been studied in several 
stereological applications (see, for instance 0). The principle used to obtain the Euler 
number of an n-dimensional domain in M" (n = 2 or 3), is based on what happens in 
an [n — l)-dimensional plane that sweeps through the domain. For domains in M?, this 
principle is based on Corollary p.4| . Now, we will extend this principle to domains in a 
surface. 

Let £) C S* C be a domain with boundary in an orientable smooth surface 5* in 
and let m G S*^. When A varies on M, the different planes tt^^a can be considered as a 
"sweeping" plane in M.^, and Corollary p.3| can be expressed as: 

xiD) = ih-B,) + ^{h-B^), 

where h, Bi denote the number of islands and bridges, respectively, observed in the 
level curves tt^ a fl S (see Figure 2) and 12,-825 similarly to [0], denote the number of 
islands and bridges observed in the "sweeping" plane 7r„ a, which contribute to the sum 
Exes(h„) Ind^/i^) (see Figure 2). 
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Figure 2 
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